Type I interferons induce an epigenetically distinct memory B cell subset in chronic viral infection.

Memory B cells (MBCs) are key providers of long-lived immunity against infectious disease, yet in chronic viral infection, they do not produce effective protection. How chronic viral infection disrupts MBC development and whether such changes are reversible remain unknown. Through single-cell (sc)ATAC-seq and scRNA-seq during acute versus chronic lymphocytic choriomeningitis viral infection, we identified a memory subset enriched for interferon (IFN)-stimulated genes (ISGs) during chronic infection that was distinct from the T-bet+ subset normally associated with chronic infection. Blockade of IFNAR-1 early in infection transformed the chromatin landscape of chronic MBCs, decreasing accessibility at ISG-inducing transcription factor binding motifs and inducing phenotypic changes in the dominating MBC subset, with a decrease in the ISG subset and an increase in CD11c+CD80+ cells. However, timing was critical, with MBCs resistant to intervention at 4 weeks post-infection. Together, our research identifies a key mechanism to instruct MBC identity during viral infection.

FixNCut: single-cell genomics through reversible tissue fixation and dissociation.

The use of single-cell technologies for clinical applications requires disconnecting sampling from downstream processing steps. Early sample preservation can further increase robustness and reproducibility by avoiding artifacts introduced during specimen handling. We present FixNCut, a methodology for the reversible fixation of tissue followed by dissociation that overcomes current limitations. We applied FixNCut to human and mouse tissues to demonstrate the preservation of RNA integrity, sequencing library complexity, and cellular composition, while diminishing stress-related artifacts. Besides single-cell RNA sequencing, FixNCut is compatible with multiple single-cell and spatial technologies, making it a versatile tool for robust and flexible study designs.

Novel Approach for Pancreas Transcriptomics Reveals the Cellular Landscape in Homeostasis and Acute Pancreatitis.

BACKGROUND & AIMS: Acinar cells produce digestive enzymes that impede transcriptomic characterization of the exocrine pancreas. Thus, single-cell RNA-sequencing studies of the pancreas underrepresent acinar cells relative to histological expectations, and a robust approach to capture pancreatic cell responses in disease states is needed. We sought to innovate a method that overcomes these challenges to accelerate study of the pancreas in health and disease. METHODS: We leverage FixNCut, a single-cell RNA-sequencing approach in which tissue is reversibly fixed with dithiobis(succinimidyl propionate) before dissociation and single-cell preparation. We apply FixNCut to an established mouse model of acute pancreatitis, validate findings using GeoMx whole transcriptome atlas profiling, and integrate our data with prior studies to compare our method in both mouse and human pancreas datasets. RESULTS: FixNCut achieves unprecedented definition of challenging pancreatic cells, including acinar and immune populations in homeostasis and acute pancreatitis, and identifies changes in all major cell types during injury and recovery. We define the acinar transcriptome during homeostasis and acinar-to-ductal metaplasia and establish a unique gene set to measure deviation from normal acinar identity. We characterize pancreatic immune cells, and analysis of T-cell subsets reveals a polarization of the homeostatic pancreas toward type-2 immunity. We report immune responses during acute pancreatitis and recovery, including early neutrophil infiltration, expansion of dendritic cell subsets, and a substantial shift in the transcriptome of macrophages due to both resident macrophage activation and monocyte infiltration. CONCLUSIONS: FixNCut preserves pancreatic transcriptomes to uncover novel cell states during homeostasis and following pancreatitis, establishing a broadly applicable approach and reference atlas for study of pancreas biology and disease.

HAMSAB diet ameliorates dysfunctional signaling in pancreatic islets in autoimmune diabetes.

An altered gut microbiota is associated with type 1 diabetes (T1D), affecting the production of short-chain fatty acids (SCFA) and glucose homeostasis. We previously demonstrated that enhancing serum acetate and butyrate using a dietary supplement (HAMSAB) improved glycemia in non-obese diabetic (NOD) mice and patients with established T1D. The effects of SCFA on immune-infiltrated islet cells remain to be clarified. Here, we performed single-cell RNA sequencing on islet cells from NOD mice fed an HAMSAB or control diet. HAMSAB induced a regulatory gene expression profile in pancreas-infiltrated immune cells. Moreover, HAMSAB maintained the expression of ß-cell functional genes and decreased cellular stress. HAMSAB-fed mice showed preserved pancreatic endocrine cell identity, evaluated by decreased numbers of poly-hormonal cells. Finally, SCFA increased insulin levels in human ß-like cells and improved transplantation outcome in NOD/SCID mice. Our findings support the use of metabolite-based diet as attractive approach to improve glucose control in T1D.

A machine learning one-class logistic regression model to predict stemness for single cell transcriptomics and spatial omics.

Cell annotation is a crucial methodological component to interpreting single cell and spatial omics data. These approaches were developed for single cell analysis but are often biased, manually curated and yet unproven in spatial omics. Here we apply a stemness model for assessing oncogenic states to single cell and spatial omic cancer datasets. This one-class logistic regression machine learning algorithm is used to extract transcriptomic features from non-transformed stem cells to identify dedifferentiated cell states in tumors. We found this method identifies single cell states in metastatic tumor cell populations without the requirement of cell annotation. This machine learning model identified stem-like cell populations not identified in single cell or spatial transcriptomic analysis using existing methods. For the first time, we demonstrate the application of a ML tool across five emerging spatial transcriptomic and proteomic technologies to identify oncogenic stem-like cell types in the tumor microenvironment.

Ligand-dependent hedgehog signaling maintains an undifferentiated, malignant osteosarcoma phenotype.

TP53 and RB1 loss-of-function mutations are common in osteosarcoma. During development, combined loss of TP53 and RB1 function leads to downregulation of autophagy and the aberrant formation of primary cilia, cellular organelles essential for the transmission of canonical Hedgehog (Hh) signaling. Excess cilia formation then leads to hypersensitivity to Hedgehog (Hh) ligand signaling. In mouse and human models, we now show that osteosarcomas with mutations in TP53 and RB1 exhibit enhanced ligand-dependent Hh pathway activation through Smoothened (SMO), a transmembrane signaling molecule required for activation of the canonical Hh pathway. This dependence is mediated by hypersensitivity to Hh ligand and is accompanied by impaired autophagy and increased primary cilia formation and expression of Hh ligand in vivo. Using a conditional genetic mouse model of Trp53 and Rb1 inactivation in osteoblast progenitors, we further show that deletion of Smo converts the highly malignant osteosarcoma phenotype to benign, well differentiated bone tumors. Conversely, conditional overexpression of SHH ligand, or a gain-of-function SMO mutant in committed osteoblast progenitors during development blocks terminal bone differentiation. Finally, we demonstrate that the SMO antagonist sonidegib (LDE225) induces growth arrest and terminal differentiation in vivo in osteosarcomas that express primary cilia and Hh ligand combined with mutations in TP53. These results provide a mechanistic framework for aberrant Hh signaling in osteosarcoma based on defining mutations in the tumor suppressor, TP53.

Fancm has dual roles in the limiting of meiotic crossovers and germ cell maintenance in mammals.

Meiotic crossovers are required for accurate chromosome segregation and producing new allelic combinations. Meiotic crossover numbers are tightly regulated within a narrow range, despite an excess of initiating DNA double-strand breaks. Here, we reveal the tumor suppressor FANCM as a meiotic anti-crossover factor in mammals. We use unique large-scale crossover analyses with both single-gamete sequencing and pedigree-based bulk-sequencing datasets to identify a genome-wide increase in crossover frequencies in Fancm-deficient mice. Gametogenesis is heavily perturbed in Fancm loss-of-function mice, which is consistent with the reproductive defects reported in humans with biallelic FANCM mutations. A portion of the gametogenesis defects can be attributed to the cGAS-STING pathway after birth. Despite the gametogenesis phenotypes in Fancm mutants, both sexes are capable of producing offspring. We propose that the anti-crossover function and role in gametogenesis of Fancm are separable and will inform diagnostic pathways for human genomic instability disorders.

Systematic benchmarking of single-cell ATAC-sequencing protocols.

Single-cell assay for transposase-accessible chromatin by sequencing (scATAC-seq) has emerged as a powerful tool for dissecting regulatory landscapes and cellular heterogeneity. However, an exploration of systemic biases among scATAC-seq technologies has remained absent. In this study, we benchmark the performance of eight scATAC-seq methods across 47 experiments using human peripheral blood mononuclear cells (PBMCs) as a reference sample and develop PUMATAC, a universal preprocessing pipeline, to handle the various sequencing data formats. Our analyses reveal significant differences in sequencing library complexity and tagmentation specificity, which impact cell-type annotation, genotype demultiplexing, peak calling, differential region accessibility and transcription factor motif enrichment. Our findings underscore the importance of sample extraction, method selection, data processing and total cost of experiments, offering valuable guidance for future research. Finally, our data and analysis pipeline encompasses 169,000 PBMC scATAC-seq profiles and a best practices code repository for scATAC-seq data analysis, which are freely available to extend this benchmarking effort to future protocols.

Human prefrontal cortex gene regulatory dynamics from gestation to adulthood at single-cell resolution.

Human brain development is underpinned by cellular and molecular reconfigurations continuing into the third decade of life. To reveal cell dynamics orchestrating neural maturation, we profiled human prefrontal cortex gene expression and chromatin accessibility at single-cell resolution from gestation to adulthood. Integrative analyses define the dynamic trajectories of each cell type, revealing major gene expression reconfiguration at the prenatal-to-postnatal transition in all cell types followed by continuous reconfiguration into adulthood and identifying regulatory networks guiding cellular developmental programs, states, and functions. We uncover links between expression dynamics and developmental milestones, characterize the diverse timing of when cells acquire adult-like states, and identify molecular convergence from distinct developmental origins. We further reveal cellular dynamics and their regulators implicated in neurological disorders. Finally, using this reference, we benchmark cell identities and maturation states in organoid models. Together, this captures the dynamic regulatory landscape of human cortical development.

Paraspeckle subnuclear bodies depend on dynamic heterodimerisation of DBHS RNA-binding proteins via their structured domains.

RNA-binding proteins of the DBHS (Drosophila Behavior Human Splicing) family, NONO, SFPQ, and PSPC1 have numerous roles in genome stability and transcriptional and posttranscriptional regulation. Critical to DBHS activity is their recruitment to distinct subnuclear locations, for example, paraspeckle condensates, where DBHS proteins bind to the long noncoding RNA NEAT1 in the first essential step in paraspeckle formation. To carry out their diverse roles, DBHS proteins form homodimers and heterodimers, but how this dimerization influences DBHS localization and function is unknown. Here, we present an inducible GFP-NONO stable cell line and use it for live-cell 3D-structured illumination microscopy, revealing paraspeckles with dynamic, twisted elongated structures. Using siRNA knockdowns, we show these labeled paraspeckles consist of GFP-NONO/endogenous SFPQ dimers and that GFP-NONO localization to paraspeckles depends on endogenous SFPQ. Using purified proteins, we confirm that partner swapping between NONO and SFPQ occurs readily in vitro. Crystallographic analysis of the NONO-SFPQ heterodimer reveals conformational differences to the other DBHS dimer structures, which may contribute to partner preference, RNA specificity, and subnuclear localization. Thus overall, our study suggests heterodimer partner availability is crucial for NONO subnuclear distribution and helps explain the complexity of both DBHS protein and paraspeckle dynamics through imaging and structural approaches.

Single-nuclei and bulk-tissue gene-expression analysis of pheochromocytoma and paraganglioma links disease subtypes with tumor microenvironment.

Pheochromocytomas (PC) and paragangliomas (PG) are rare neuroendocrine tumors associated with autonomic nerves. Here we use single-nuclei RNA-seq and bulk-tissue gene-expression data to characterize the cellular composition of PCPG and normal adrenal tissues, refine tumor gene-expression subtypes and make clinical and genotypic associations. We confirm seven PCPG gene-expression subtypes with significant genotype and clinical associations. Tumors with mutations in VHL, SDH-encoding genes (SDHx) or MAML3-fusions are characterized by hypoxia-inducible factor signaling and neoangiogenesis. PCPG have few infiltrating lymphocytes but abundant macrophages. While neoplastic cells transcriptionally resemble mature chromaffin cells, early chromaffin and neuroblast markers are also features of some PCPG subtypes. The gene-expression profile of metastatic SDHx-related PCPG indicates these tumors have elevated cellular proliferation and a lower number of non-neoplastic Schwann-cell-like cells, while GPR139 is a potential theranostic target. Our findings therefore clarify the diverse transcriptional programs and cellular composition of PCPG and identify biomarkers of potential clinical significance.

SSNIP-seq: A simple and rapid method for isolation of single-sperm nucleic acid for high-throughput sequencing.

We developed a simple and reliable method for the isolation of haploid nuclei from fresh and frozen testes. The described protocol uses readily available reagents in combination with flow cytometry to separate haploid and diploid nuclei. The protocol can be completed within 1 hour and the resulting individual haploid nuclei have intact morphology. The isolated nuclei are suitable for library preparation for high-throughput DNA and RNA sequencing using bulk or single nuclei. The protocol was optimised with mouse testes and we anticipate that it can be applied for the isolation of mature sperm from other mammals including humans.

Inhibition of mutant IDH1 promotes cycling of acute myeloid leukemia stem cells.

Approximately 20% of acute myeloid leukemia (AML) patients carry mutations in IDH1 or IDH2 that result in over-production of the oncometabolite D-2-hydroxyglutarate (2-HG). Small molecule inhibitors that block 2-HG synthesis can induce complete morphological remission; however, almost all patients eventually acquire drug resistance and relapse. Using a multi-allelic mouse model of IDH1-mutant AML, we demonstrate that the clinical IDH1 inhibitor AG-120 (ivosidenib) exerts cell-type-dependent effects on leukemic cells, promoting delayed disease regression. Although single-agent AG-120 treatment does not fully eradicate the disease, it increases cycling of rare leukemia stem cells and triggers transcriptional upregulation of the pyrimidine salvage pathway. Accordingly, AG-120 sensitizes IDH1-mutant AML to azacitidine, with the combination of AG-120 and azacitidine showing vastly improved efficacy in vivo. Our data highlight the impact of non-genetic heterogeneity on treatment response and provide a mechanistic rationale for the observed combinatorial effect of AG-120 and azacitidine in patients.

Isolating Nuclei From Frozen Human Heart Tissue for Single-Nucleus RNA Sequencing.

Heart disease is the leading cause of global morbidity and mortality. This is in part because, despite an abundance of animal and in vitro models, it has been a challenge to date to study human heart tissue with sufficient depth and resolution to develop disease-modifying therapies for common cardiac conditions. Single-nucleus RNA sequencing (snRNA-seq) has emerged as a powerful tool capable of analyzing cellular function and signaling in health and disease, and has already contributed to significant advances in areas such as oncology and hematology. Employing snRNA-seq technology on flash-frozen human tissue has the potential to unlock novel disease mechanisms and pathways in any organ. Studying the human heart using snRNA-seq is a key priority for the field of cardiovascular sciences; however, progress to date has been slowed by numerous barriers. One key challenge is the fact that the human heart is very resistant to shearing and stress, making tissue dissociation and nuclear isolation difficult. Here, we describe a tissue dissociation method allowing the efficient and cost-effective isolation of high-quality nuclei from flash-frozen human heart tissue collected in surgical operating rooms. Our protocol addresses the challenge of nuclear isolation from human hearts, enables snRNA-seq of the human heart, and paves the way for an improved understanding of the human heart in health and disease. Ultimately, this will be key to uncovering signaling pathways and networks amenable to therapeutic intervention and the development of novel biomarkers and disease-modifying therapies. © 2022 Wiley Periodicals LLC. Basic Protocol: Human heart tissue dissociation and nuclear isolation for snRNA-seq.

ALTEN: A High-Fidelity Primary Tissue-Engineering Platform to Assess Cellular Responses Ex Vivo.

To fully investigate cellular responses to stimuli and perturbations within tissues, it is essential to replicate the complex molecular interactions within the local microenvironment of cellular niches. Here, the authors introduce Alginate-based tissue engineering (ALTEN), a biomimetic tissue platform that allows ex vivo analysis of explanted tissue biopsies. This method preserves the original characteristics of the source tissue's cellular milieu, allowing multiple and diverse cell types to be maintained over an extended period of time. As a result, ALTEN enables rapid and faithful characterization of perturbations across specific cell types within a tissue. Importantly, using single-cell genomics, this approach provides integrated cellular responses at the resolution of individual cells. ALTEN is a powerful tool for the analysis of cellular responses upon exposure to cytotoxic agents and immunomodulators. Additionally, ALTEN's scalability using automated microfluidic devices for tissue encapsulation and subsequent transport, to enable centralized high-throughput analysis of samples gathered by large-scale multicenter studies, is shown.

Epigenetic Activation of Plasmacytoid DCs Drives IFNAR-Dependent Therapeutic Differentiation of AML.

Pharmacologic inhibition of epigenetic enzymes can have therapeutic benefit against hematologic malignancies. In addition to affecting tumor cell growth and proliferation, these epigenetic agents may induce antitumor immunity. Here, we discovered a novel immunoregulatory mechanism through inhibition of histone deacetylases (HDAC). In models of acute myeloid leukemia (AML), leukemia cell differentiation and therapeutic benefit mediated by the HDAC inhibitor (HDACi) panobinostat required activation of the type I interferon (IFN) pathway. Plasmacytoid dendritic cells (pDC) produced type I IFN after panobinostat treatment, through transcriptional activation of IFN genes concomitant with increased H3K27 acetylation at these loci. Depletion of pDCs abrogated panobinostat-mediated induction of type I IFN signaling in leukemia cells and impaired therapeutic efficacy, whereas combined treatment with panobinostat and IFNα improved outcomes in preclinical models. These discoveries offer a new therapeutic approach for AML and demonstrate that epigenetic rewiring of pDCs enhances antitumor immunity, opening the possibility of exploiting this approach for immunotherapies. SIGNIFICANCE: We demonstrate that HDACis induce terminal differentiation of AML through epigenetic remodeling of pDCs, resulting in production of type I IFN that is important for the therapeutic effects of HDACis. The study demonstrates the important functional interplay between the immune system and leukemias in response to HDAC inhibition. This article is highlighted in the In This Issue feature, p. 1397.

Same-Cell Co-Occurrence of RAS Hotspot and BRAF V600E Mutations in Treatment-Naive Colorectal Cancer.

PURPOSE: Mitogen-activated protein kinase pathway-activating mutations occur in the majority of colorectal cancer (CRC) cases and show mutual exclusivity. We identified 47 epidermal growth factor receptor/BRAF inhibitor-naive CRC patients with dual RAS hotspot/BRAF V600E mutations (CRC-DD) from a cohort of 4,561 CRC patients with clinical next-generation sequencing results. We aimed to define the molecular phenotypes of the CRC-DD and to test if the dual RAS hotspot/BRAF V600E mutations coexist within the same cell. MATERIALS AND METHODS: We developed a single-cell genotyping method with a mutation detection rate of 96.3% and a genotype prediction accuracy of 92.1%. Mutations in the CRC-DD cohort were analyzed for clonality, allelic imbalance, copy number, and overall survival. RESULTS: Application of single-cell genotyping to four CRC-DD revealed the co-occurrence of both mutations in the following percentages of cells per case: NRAS G13D/KRAS G12C, 95%; KRAS G12D/NRAS G12V, 48%; BRAF V600E/KRAS G12D, 44%; and KRAS G12D/NRAS G13V, 14%, respectively. Allelic imbalance favoring the oncogenic allele was less frequent in CRC-DD (24 of 76, 31.5%, somatic mutations) compared with a curated cohort of CRC with a single-driver mutation (CRC-SD; 119 of 232 mutations, 51.3%; P = .013). Microsatellite instability-high status was enriched in CRC-DD compared with CRC-SD (23% v 11.4%, P = .028). Of the seven CRC-DD cases with multiregional sequencing, five retained both driver mutations throughout all sequenced tumor sites. Both CRC-DD cases with discordant multiregional sequencing were microsatellite instability-high. CONCLUSION: Our findings indicate that dual-driver mutations occur in a rare subset of CRC, often within the same tumor cells and across multiple tumor sites. Their presence and a lower rate of allelic imbalance may be related to dose-dependent signaling within the mitogen-activated protein kinase pathway.

Non-genetic determinants of malignant clonal fitness at single-cell resolution.

All cancers emerge after a period of clonal selection and subsequent clonal expansion. Although the evolutionary principles imparted by genetic intratumour heterogeneity are becoming increasingly clear1, little is known about the non-genetic mechanisms that contribute to intratumour heterogeneity and malignant clonal fitness2. Here, using single-cell profiling and lineage tracing (SPLINTR)-an expressed barcoding strategy-we trace isogenic clones in three clinically relevant mouse models of acute myeloid leukaemia. We find that malignant clonal dominance is a cell-intrinsic and heritable property that is facilitated by the repression of antigen presentation and increased expression of the secretory leukocyte peptidase inhibitor gene (Slpi), which we genetically validate as a regulator of acute myeloid leukaemia. Increased transcriptional heterogeneity is a feature that enables clonal fitness in diverse tissues and immune microenvironments and in the context of clonal competition between genetically distinct clones. Similar to haematopoietic stem cells3, leukaemia stem cells (LSCs) display heritable clone-intrinsic properties of high, and low clonal output that contribute to the overall tumour mass. We demonstrate that LSC clonal output dictates sensitivity to chemotherapy and, although high- and low-output clones adapt differently to therapeutic pressure, they coordinately emerge from minimal residual disease with increased expression of the LSC program. Together, these data provide fundamental insights into the non-genetic transcriptional processes that underpin malignant clonal fitness and may inform future therapeutic strategies.